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A short survey of the semiclassical periodic orbit theory, initiated by M. Gutzwiller and general- 
ized by many other authors, is given. Via so-called semiclassical trace formmulae, gross-shell effects 
in bound fermion systems can be interpreted in terms of a few periodic orbits of the corresponding 
classical systems. In integrable systems, these are usually the shortest members of the most degen- 
erate families or orbits, but in some systems also less degenerate orbits can determine the gross-shell 
structure. Applications to nuclei, metal clusters, semiconductor nanostructures, and trapped dilute 
atom gases are discussed. 
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I. INTRODUCTION 

The Nilsson model Q gave rise to my first steps 
towards scientific research. My teacher in theoretical 
physics was Kurt Alder, who had been a member of 
the theoretical 'Coulomb excitation crew' at the Niels 
Bohr Institute during the years when collective nuclear 
motion was explored and single-particle motion in de- 
formed nuclei was studied. When I asked Alder about 
a possible subject for my diploma work, he gave me 
a copy of Sven Gosta's famous paper and asked me 
to write a program for computing the deformed single- 
particle wavefunctions; he had become weary of interpo- 
lating their expansion coefficients in Nilsson's tables. I 
set about to diagonalize the Nilsson Hamiltonian using, 
however, not the spherical basis but that of the deformed 
harmonic oscillator, i.e., the asymptotic Nilsson states. 
The results were bound to be the same as those of Nils- 
son and were therefore not published j3l - but through 
this exercise I became initiated into the shell structure of 
nuclei. I felt very privileged when I later came to know 
the great human being behind the famous model. 

In this paper I want to review the semiclassical de- 
scription of shell-effects in finite fermion systems using 
the periodic orbit theory. After a brief reminder about 
trace formulae, I will discuss some of their applications 
to systems in four different branches of physics. Since all 
results have been published elsewhere, I will not repro- 
duce here any figures, but just discuss in words some of 
the most important results and conclusions. 

II. PERIODIC ORBIT THEORY AND 
GROSS-SHELL EFFECTS 

A. Semiclassical trace formulae 

The periodic orbit theory (POT) was initiated by M. 
Gutzwiller in a series of publications culminating in his 
seminal paper in 1971 that contains the semiclassical 
trace formula |^]. It relates the quantum spectrum {Ei} 
(which we here assume to be discrete, although the in- 
clusion of a continuum is possible) of a hermitian Hamil- 



tonian H to the periodic orbits of the corresponding 
classical Hamiltonian H{q, p). The quantum-mechanical 
level density, defined as the sum of Dirac delta functions 
peaked at the levels Ei, can be decomposed into a smooth 
part g{E) and an oscillating part 5g{E): 

g{E) = Y,5{E~E,)^g{E) + 5g{E). (1) 

i 

The smooth part contains by definition the average level 
density which usually is a monotonously increasing func- 
tion of E and can be obtained in the extended Thomas- 
Fermi (ETF) model (see, e.g., ^;5], Ch. 4). The oscillating 
part can be expressed by the semiclassical trace formula 

5gsc{E) ~ \Sj,o{E)/n - apo7r/2] . (2) 

po 

The sum is over all periodic orbits (po) of the classical 
system, Spo{E) — f p ■ dq are their action integrals, the 
amplitudes Apo{E) depend on their stabilities and de- 
generacies, and tJpo are called the Maslov indices. The 
sum in ^ is an asymptotic one, correct to leading order 
in 1/h, and in non-integrable systems it is hampered by 
convergence problems 6:j. For isolated orbits, Gutzwiller 
expressed Q the amplitudes Apo{E) in terms of their 
periods and stability matrices. 

The trace formula Q was later generalized to billiard 
syster ns iTj and to systems with continuous symmetries 
H, 0, including integrable systems. A relativistic 
trace formula for spin 1/2 particles was derived in |l l| , 
and a nonrelativistic trace formula for particles with arbi- 
trary spin s in |0 . In all cases, the trace formula has the 
same general form 0, but the amplitudes Apo{E) take 
different forms. For isolated orbits, their h dependence is 
given by a factor h~ ^ , while for orbits appearing in /-fold 
degenerate families, the amplitudes go like . 

In integrable and mixed-dynamical systems, periodic 
orbits can change their stability under the variation of a 
control parameter (e.g., the energy E, a potential param- 
eter, or an additional external field) and thereby undergo 
bifurcations. In such situations, the amplitudes Apo di- 
verge at the bifurcation points. The same happens also 
in limits where continuous symmetries are broken (or re- 
stored), since hereby the h dependence of the Apo changes 
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discontinuously. The remedy to remove these (unphys- 
ical!) divergences is to go beyond the stationary-phase 
approximation for the integration(s) used in the deriva- 
tion of the semiclassical trace formula. This has, besides 
0, H, J been devefoped most systematically in fi^ for 
symmetry breaking and bifurcations, and in |l4| for sym- 
metry breaking in weakly perturbed integrable systems, 
leading in all cases to local uniform approximations with 
finite amplitudes Apo- Global uniform approximations 
which yield finite amplitudes at symmetry-breaking and 
bifurcation points, and far from them go over into the 
standard (extended) Gutzwiller trace formula, were de- 
veloped for the breaking of U(l) symmetry in for 
some cases of U(2) and S0(3) symmetry breaking in [Tfil |. 
for the symmetry breaking U(3) -> SO (3) in Tt'I, and for 
various types of bifurcations in .ISj . (Details and further 
references may be found in 0, Ch. 6.3.) 

For interacting finite fermion systems described in the 
mean- field approximation (i.e., in Hartree-Fock or den- 
sity functional theory), one can also obtain semiclassical 
trace formulae for the oscillating parts of the total bind- 
ing energy Eh and the particle number N. Hereby one 
writes, similarly to Eb = Eb + SE and N = N + SN. 
The average quantities Eb and N are taken from the ETF 
model, and for the oscillating parts one finds 0, @ 



SEsc^ V^po(A) 
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both to be evaluated at the Fermi energy X{N) for a given 
number of particles A'^. The periodic orbits are ideally 
those of the classical counterpart of the self-consistent 
mean field, which for practical purposes often is taken as 
a shell-model type potential. 



B. Coarse-graining and finite temperatures 

Our present emphasis in the use of POT is not the 
full quantization of the spectra of finite fermion systems, 
but on the semiclassical description of their gross-shell 
structure. For this purpose we coarse-grain the quantum 
spectrum by a convolution of the level density Q with a 
normalized Gaussian of width 7: 
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The coarse graining of the trace formula ^ gives, using 
the stationary-phase approximation for the convolution 
integral, an extra exponential factor in the trace formula: 
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The same exponential factor appears in the trace formu- 
lae ^ for SEsc{X,j) and SNsdX,^)- It suppresses the 
contributions from orbits with larger periods Tpo. A simi- 
lar suppression of longer orbits and the overall amplitude 
of the shell effects occurs at finite temperatures. E.g. in a 
grand-canonical system at temperature T, the oscillating 
part of the free Helmholtz energy has the trace formula 
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Here Tpo — kBTTrTpo{X)/h and ks is the Boltzmann con- 
tant. In both situations the gross-shell effects are domi- 
nated by the shortest periodic orbits of the system ^ . 

In mixed-dynamical and integrable systems, orbits 
with different degrees / of degeneracy can coexist. The 
gross-shell structure then results from a competition be- 
tween the periods and the degeneracies of the shortest 
orbits Whereas their coarse-grained semiclassical 

amplitudes decrease with growing period Tpo, they in- 
crease with growing degeneracy / due to their depen- 
dence oc already mentioned above. 

In arbitrary spherical three-dimensional systems, the 
most degenerate orbits undergo both radial and angular 
oscillations with rational frequency ratios (cf. [fll Il7ll2(l | ) : 



e N. 
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(For physical reasons, only pairs of integers n, m with 
equal signs are allowed; orbits with negative n, m cor- 
respond to the time-reversed of the orbits with positive 
n,m). The orientiations of these 'rational tori' can be 
rotated about three Euler angles without changing their 
shapes, periods or actions; therefore they appear in three- 
fold degenerate families (/ = 3). Their existence for ar- 
bitrary ratios n:m depends, however, on the form of the 
radial potential V{r) (cf. [l7ll21| ) and, in general, on the 
energy E. The special orbits with angular momentum 
L — and L — Lmax , corresponding to librating 'diame- 
ter' and rotating 'circle' orbits, respetively, form families 
with only two- fold degeneracy (/ = 2), since one of the 
three Euler rotations does not change their orientations. 
The Coulomb potential V(r) = —a/r has an extra dy- 
namical symmetry, leading to 0(4); here all orbits have 
n:m — 1:1. Spherical harmonic oscillators V{r) = ar"^ 
also have an extra dynamical symmetry, leading in three 
dimensions to U(3); here all orbits have n : m = 2:1. 
In each of these two special potentials, all periodic or- 
bits (including the diameters and circles) form one fam- 
ily with degeneracy / = 4. Quantum-mechanically, the 
dynamical symmetries reflect themselves in an accidental 
extra degeneracy of the eigenvalue spectrum {Ei}. The 
semiclassical trace formulae for Sg{E) of these systems, 
added to their ETF expressions for g (E) , reproduce the 
exact quantum-mechanical level densities according to 
(QJ [see Ji], Eqs. (3.144) and (3.69) for their explicit an- 
alytical expressions]. 
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III. POT FOR FINITE FERMION SYSTEMS 
A. Nuclei 

1. Ground-state deformations. Strutinsky et al. [^l22|| 
were the first to extend the POT to systems with contin- 
uous symmetries and to apply it to the study of gross- 
shell effects in nuclei. In |23 they studied the periodic 
orbits in a spheroidal cavity with axis ratio 77 as a model 
for the mean field of a deformed nucleus. They plot- 
ted the shell-correction energy SE, obtained quantum- 
mechanically using Strutinsky's shell-correction method 
psf with the spectra of realistic deformed Woods-Saxon 
potentials |3| with the same spheroidal deformations (in- 
cluding spin-orbit interaction), versus particle number N 
and deformation ry. The slopes of the valleys in these de- 
formation energy surfaces SE(N, rj) corresponding to the 
ground-state deformations could then be correctly repro- 
duced by the condition that the actions Spo of the short- 
est and most degenerate periodic orbits in the spheroidal 
cavity (with f = 2) j25| be constant. Contributions of 
orbits with / = 1 were negligible. (The Fermi energies 
corresponding to the spherical magic numbers had to be 
adjusted, as no spin-orbit interaction was included in the 
cavity model.) Although this was only a qualitative re- 
sult, it proved the correctness of the concept to interpret 
quantum-mechanical gross-shell effects semiclassically in 
terms of short periodic orbits of the corresponding clas- 
sical system. The use of a cavity with infinitely steep 
walls for the mean field hereby justifies itself through the 
short range and the saturating property of the effective 
nucleon-nucleon interaction, leading to steep walls of the 
self-consistent Hartree-Fock potentials or their approxi- 
mations by Woods-Saxon type shell-model potentials. 

2. Left-right asymmetry of fission barriers. A promi- 
nent manifestation of shell effects is the 'double-humped' 
fission barrier of nuclei in the actinide region . One 
particular aspect is that of the onset of a left-right asym- 
metry of the fissioning nuclear shapes which eventually 
leads to the asymmetric mass distributions of the fission 
fragments. Since Sven Gosta Nilsson and his group, and 
other scientists in Lund, were much involved in the study 
of this shell effect, I may dwell a little on its history. 
The mixing of pairs of single-particle states with oppo- 
site parities in a spheroidal harmonic-oscillator potential 
was studied earlier as a possible mechanism lead- 
ing to 'pear-shaped' nuclei. In 1962, S. A. E. Johansson 
[23 took this question up and investigated the possibil- 
ity of octupole-deformed fission barriers. Since Strutin- 
sky's shell-correction method 23] did not yet exist at 
that time, no realistic fission barriers could be obtained 
qantum-mcchanically with the Nilsson model. Johans- 
son showed, however, that at the typical deformations of 
the actinide fission barriers predicted by the liquid-drop 
model (LDM) |29l| . th e mixing of single-particle orbits of 
the type used in [23| leads to an instability of the bar- 
rier against octopole shapes. Using the shell-correction 
method with the Nilsson model, P. MoUer and S. G. Nils- 



son [s^l obtained in 1970 the instability of the outer fis- 
sion barrier against a suitable mixture of £3 and £5 defor- 
mations. Thus, the onset of the fission mass asymmetry 
was clearly a quantum-mechanical shell effect that could 
not be explained by the classical LDM model. In a de- 
tailed microscopical study, C. Gustafsson, P. MoUer, and 
S. G. Nilsson '3ll| showed a year later that those pairs 
of single-particle states, which are most sensitive to the 
left-right asymmetric shapes and hence responsible for 
their onset, have their wavefunction nodes and extrema 
on parallel planes at and near the waist-line of the fis- 
sioning nucleus perpendicular to its symmetry axis. 
30 years later, in a Lund-Regensburg-Dresden collabo- 
ration this effect was studied semiclassically. The 
POT had been used in the same collaboration for 
cavities with the {c,h,a) shapes of 24], for which the 
mass asymmetry of the outer fission barrier of actinide 
nuclei had also been obtained quantum-mechanically in 
[33 |. The shortest periodic orbits here are families with 
/ = 1, having the axial U(l) symmetry; they are sim- 
ply the diagonal, triangular and square-shaped orbits in 
the circular planes perpendicular to the nuclear symme- 
try axis. The semiclassical trace formula (with a uniform 
approximation for the bifurcation occurring at the onset 
of the neck, where the orbits in the central plane be- 
come unstable and give birth to two new parallel planes 
with stable orbits) was shown [3^ to yield realistic de- 
formation encry surfaces SEsc{c,h,a) in the region of 
the outer fission barrier, predicting its instability against 
the asymmetry parameter a in good agreement with the 
old quantum-mechanical results [s^. (Again, the spin- 
orbit interaction was omitted; the Fermi energy as the 
only parameter was adjusted to yield the fission iso- 
mer minimum at the correct deformation.) Similarly to 
[2^ . the valley of steepest descent through the deforma- 
tion energy surface, leading over a left-right asymmet- 
ric outer saddle, is obtained by the stationary condition 
SSpo = for the shortest orbits. In 32], it was shown 
that these orbits are situated in a very small regular is- 
land of a dominantly chaotic phase space. An approxi- 
mate Einstein-Brillouin-Keller (EBK) quantization of the 
linearized classical motion in these regular islands repro- 
duced rather precisely the quantum-mechanical energies 
of those diabatic single-particle states with opposite par- 
ity which are most sensitive to the a deformations and 
hence quantum-mechanically responsible for the a insta- 
bility of the outer barrier. Furthermore, their wavefunc- 
tions were found to have their nodes and extrema pre- 
cisely in the planes near the nuclear waist-line that con- 
tain the shortest periodic orbits responsible semiclassi- 
cally for the asymmetry effect. 

This application of the POT represents an interesting ex- 
ample for the classical-to-quantum correspondence of the 
interplay between chaos and order: a tiny regular island 
in an almost chaotic phase space causes a quantum shell 
effect in an interacting many-body system with observ- 
able consequences in the form of the mass asymmetry of 
the fission fragments. 
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B. Metal clusters 

Metal clusters are interesting finite fermion systems 
which allow one to study the transition from atoms over 
molecules towards condensed matter 35j . In the simplest 
theoretical description, the so-called 'jellium model', the 
ions are replaced by a structureless but deformable posi- 
tive background and the systems of N interacting valence 
electrons in the external jellium potential are studied |3(ij . 
In neutral clusters, the eletrostatic long-range forces can- 
cel and the valence electrons are only bound by the short- 
ranged exchange and correlation effects. Neutral metal 
clusters therefore have much in common with nuclei. One 
difference to nuclei is that there is no measurable spin- 
orbit interaction in most metal clusters. The magic num- 
bers Ni of the smallest spherically stable clusters corre- 
spond to those of a harmonic oscillator {Ni — 2, 8, 20, 40) 
|85| ■ For the analysis of early experimental abundance 
spectra of small sodium clusters, the Nilsson model with- 
out spin-orbit term was successfully employed |^ to in- 
terpret the regions between the spherical shell closures 
in terms of prolate and oblate deformations. In a self- 
consistent mean- field description [ssL Is^ . the average 
potential of clusters with N ^ 80 valence electrons has 
steep walls like heavy nuclei, and therefore cavity models 
provide again a good approximation. 

One early result of POT was the observation by Balian 
and Bloch that the coarse-grained level density of a 
spherical cavity exhibits a pronounced beating pattern: a 
rapid regular oscillation, reflecting the shell structure of 
the spectrum, modulated by a slow oscillation reaching 
over some 13 - 14 shells. From the trace formula derived 
in one sees that the beat comes about by the interfer- 
ence of the shortest periodic orbits of highest degeneracy 
(/ = 3), which here are the triangle (n :m = 3:1) and 
square (4 : 1) orbits. (The diameter orbit with degener- 
acy / = 2 can be neglected.) The rapid shell oscillations 
are determined by the average length of these two orbits, 
while the period of their amplitude modulation, the so- 
called 'super-shell' oscillation, is given by the difference 
of their lengths. 

The numbers of fermions needed to reach the first 
super-shell node is of the order of 800 — 1000. Super- 
shells can therefore not be seen in nuclei. Neutral metal 
clusters, however, can be made arbitrarily large. This in- 
spired Nishioka et al. ^3 to study the super-shell struc- 
ture in Woods-Saxon potentials fitted to self-consistent 
mean fields [sl, , and to predict that it should be ob- 
servable in metal clusters. Indeed, the super-shells were 
experimentally observed, for the first time in supersonic 
beams of hot sodium clusters ji^]. Their abundance in 
an adiabatically expanding beam is dominated by their 
stability against evaporation of single atoms and exhibits 
pronounced peaks at the spherically magic numbers. The 
larger ones, Ni = 92, 138, 192, 264, . . . , are almost exactly 
those of a spherical cavity. By a suitable extraction of 
the oscillating part of the abundance s pec tra, the super- 
shell beat could clearly be exhibited |43, with its first 



node appearing around N ~ 900. 

1/3 

Plotting the cube roots of the magic numbers, N^ 
which are proportional to the r.m.s. radii of the magic 
clusters, against the shell number i, one obtains a straight 

1/3 1/3 

line with a slope s = N^^^^ ~ ^i ■ The experimental 
value of this slope is Sexp — 0.61 ± 0.01; it was confirmed 
by later experiments, also with other types of metal clus- 
ters (their different Wigner-Seitz radii do not affect the 
value of s, as long as they are single- valenced) [35ll3^. 
Around i ~ 14, there is a slight discontinuity in the 
plot s(i) before it continues again with the same slope. 
This is due to the phase change of the rapid shell oscil- 
lations when passing through the first super-shell node. 
The value of this slope predicted by the POT of 7, 4(| 
is Spot = 0.603, that of the self-consistent quantum- 
mechanical mean- field calculations is Smf = 0.61, both 
in perfect quantitative agreement with experiment. 

This provides another example of the good agreement 
of POT with both experiment and quantum mechanics. 
An easily readable account of the super-shells in metal 
clusters is given in [43| . 

Unfortunately, it is very difficult to extract spectro- 
scopic data on the electronic single-particle states in 
metallic clusters. The most direct access to their shapes 
is given by the so-called 'Mie plasmons', the collective 
dipole oscillations of the valence electrons against the 
ions 0,E3 (which are the origin of the colors in stained- 
glass windows ^3)- Similarly to the giant-dipole res- 
onances in nuclei, their splitting gives evidence of the 
average deformations of the clusters. This gave rise to a 
series of theoretical investigations of cluster deformations 
using deformed shell models (see js^). For applications 
of the POT one employs hereby most comfortably the 
spheroidal cavity model used also in @ . I refer to ^3 for 
a detained account containing also a nice application of 
the POT to the semiclassical interpretation of moments 
of inertia. The effects of weak magnetic fields on spheri- 
cal metal clusters were studied in ^3 ■ S3 > spheroidal 
cavities with the lowest multipole deformations £2, £3, 
an d ca were studied, and the perturbative trace formula 
of [g (cf. also [13) was used to predict their ground- 
state deformations. The results were in very good agree- 
ment with those of quantum-mechanical shell-correction 
calculations using the spectra of the same cavities. 

C. Semiconductor nanostructures 

Semiconductor heterostructures can be used to con- 
struct two-dimensional systems of quasi-free electrons on 
the nanometer scale. With the help of external metallic 
gates or lithography, the electrons can further be lat- 
erally confined to form so-called quantum dots, quan- 
tum channels, quantum wires, antidot superlattices, etc. 
[53, 1^- Applying a perpendicular magnetic field B, one 
can measure the magneto-resistance of such devices. Un- 
der suitable experimental circumstances, both the mean 
free path and the phase coherence length can be made 
larger than the sizes of these structures, so that quan- 
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turn interference still takes place while the dimensions 
are large enough to allow for a semiclassical description. 
Nanostructues are therefore ideal tools to study the in- 
terplay between classical and quantum mechanics. 

Weiss et al. [H^l measured the resistance of andidot su- 
perlattices and found oscillations which can be explained 
classically by the commensurability of cyclotron orbits 
with the superdot lattice: when an electron is trapped in 
a cyclotron orbit that fits around 1, 2, 4, 9, etc. antidots, 
it does not contribute to the conductance and hence a 
peak is seen in the magneto-resistance (see js^ for an 
easily readable account). In weak B fields at very low 
temperatures, some rapid B-periodic oscillations could 
be observed. They could be interpreted semiclassically 
by the interference of different trapped periodic orbits 
of comparable lenghts; the linear response of the system 
to the B field was hereby described by a semiclassical 
version of the Kubo theory yielding a trace formula for 
the conductance [H^l- These so-called 'Aharonov-Bohm 
(AB) oscillations' can be easily understood in a pertur- 
bative approach in which the effect of the magnetic 
field is taken into account to lowest order only in the 
actions Spo = § p ■ dq, while the shapes of the orbits 
are left unchanged. Under the canonical substitution 
p ^ p — (e/c)A, where A is the vector potential with 
B=:Vx A, the action of a periodic orbit changes like 

Spo Spo - (e/c) (j) A-dq = Spo - (e/c)$po , (8) 

where $po = ^ B • dFpo is the magnetic flux through the 
area Fpo surrounded by the orbit. Consequently, the 
perturbed semiclassical trace formula is modified only 
by a factor cos [e^po/hc) containing the Aharonov-Bohm 
phase which causes the _B-periodic oscillations. 

Similar AB oscillations have also been measured in the 
magneto-conductance of a circular quantum dot contain- 
ing some ^ 1200 to 2000 electrons jSH. They could be 
qualitatively well explained j56j by the perturbed level 
density Sg{E, B) of a two-dimensional circular billiard. 
This is an integrable system whose trace formula is an- 
alytically known [s^l- The conductance oscillations as 
functions of the radius of the quantum dot (regulated 
experimentally by the applied gate voltage) were well re- 
produced by the average length of the shortest orbits 



(here: diameters and triangles), while the period of the 
AB oscillations according to ((SJl was well reproduced by 
the area enclosed by the triangular orbit |5g . 

An analytical trace formula for the circle billiard in ar- 
bitrarily strong transverse magnetic fields B was given in 
|58| , and the magnetization of quantum dots was studied 
semiclassically in [59l |. 

In a mesoscopic semiconductor channel with two anti- 
dots, the magneto-conductance was measured in [g^ and 
also found to exhibit AB oscillations. These could be well 
explained using the semiclassical Kubo formula [HJ] 
with a suitably modeled two-dimensional confinement 
potential for the channel (with antidots), which repre- 
sents a mixed-dynamical system with a rather chaotic 
phase space. Plotting the maxima of the experimental 
AB oscillations versus magnetic field B and gate voltage 
(which regulates the radii of the antidots) , one obtains a 
grid of lines exhibiting some characteristic displacements 
|6Cll |. At first sight, these might be attributed to miss- 
ing flux units. Quantum-mechanical calculations [G^I re- 
produced these displacements but could not explain the 
physics behind them. In the semiclassical calculations 
the displacements could, in fact, be attributed to bifur- 
cations of some of the trapped periodic orbits [6ll l63 | . 

D. Trapped dilute atomic gases 

I mention only briefly the finite fermion systems pro- 
duced by confining diluted fermionic atom gases, e.g. 
in magneto-optic traps [B^I- Hartree-Fock calculations 
for N harmonically trapped atoms with a short-ranged 
repulsive interaction |63 yielded shell effects 6E{N) in 
their total binding energies Eb{N) which remind about 
the super-shells discussed above. More details are given 
in 66J; let me just emphasize here that the origin of the 
beating shell structure is different here from that in a 
spherical cavity 0- The self-consistent mean field can 
be modeled by a perturbed harmonic oscillator V(r) = 
_|_ gj,4_ Pqj. gm^h potentials it was shown recently 
[r^ that the gross-shell structure is dominated by the 
two-fold degenerate diameter and circle orbits, whose in- 
terference explains the super-shells found in |65j . The 
shortest three-fold degenerate orbits have frequency ra- 
tios n -.m > 7:3 and contribute only to finer details of 
the quantum spectrum at relatively high energies. 
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